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Poynting-Robertson Effect on the Lyapunov Stability of Equilibrium 
Points in the Generalized Photogravitational Chermnykh's Problem 
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Abstract The Poynting-Robertson(P-R) effect on Lyapunov stability of equilibrium points is 
being discussed in the the generalized photogravitational Chermnykh's problem when bigger 
primary is a sours of radiation and smaller primary is an oblate spheriod. We derived the 
equations of motion, obtained the equilibrium points and examined the linear stability of 
the equilibrium points for various values of parameter which have been used in the present 
problem. We have examined the effect of gravitational potential from the belt. The positions 
of the equilibrium points are different from the position in classical case. We have seen that 
due to the P-R effect all the equilibrium points are unstable in Lyapunov sense. 

Key words: equilibrium points: generalized photogravitationakChermnykh's prob- 
lem:radiation pressure: Poynting-Robertson effect 

1 INTRODUCTION 



The Chermnykh's 



Chermnvkh 



(1987). 



problem is new kind of restricte d three body problem which was first time studied by 



Papadakis and Kanav os (2007) given numerical exploration of Chermnykh's problem, 
in which the equilibrium points and zero velocity curves studied numerically also the non-linear stability 
for the triangular Lagrangian points are computed numerically for the Earth-Moon and Sun- Jupiter mass 
distribution when the angular velocity varies. The solar radiation pressure force F p is exactly apposite to 
the gravitational attraction force F g and change with the distance by the same law it is possible to consider 
that the result of action of this force will lead to reducing the effective mass of the Sun or particle. It is 
acceptable to speak about a reduced mass of the particle as the effect of reducing its mass depe nds on the 



properties of the 



(1903) 



Robertson 



^ articl e itself. The first order in — the radiation pressure force is given by[see lPovnting 



(1937)]: 



P \R 



V.RR 

cR 2 



V_ 
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Where F p = 16 ^T psc denotes the measure of the radiation pressure force, R the position vector of P with 
respect to radiation sours Sun S, V the corresponding velocity vector and c the velocity of light. In the 
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expression of F p , L is luminosity of the radiating body, while m, p and s are the mass, density and cross 
section of the particle respectively. 

The first term in equation (Q~|i expresses the radiation pressure. The second term represents the Doppler 
shift of the incident radiation and the third term is due to the absorption and subsequent re-emission of the 
incident radiation. These last two terms taken together are the Poynting-Robertson effect. The Poynting- 
Robe rtson effec t will o perate to sweep small particles of the solar system into the Sun at cosmically rapid 
rate . IChernikovl i 1 9701) discussed the position as well as the stabi l ity of the Lagrangian equilibrium points 



when radiation pressure, P-R drag force are included. iMurrayl (11994) systematicall y discussed the dy 



namic al effect of general drag in the planar circular restricted three body problem. 



Ishwar and Kushvah 



(2006) examined the linear stability of triangular equilibrium points in the generalized photogravitational 
restricted three body problem with Poynting-Robertson drag, L4 and L$ points became unstable due to P-R 
drag which is very remarkable and important, where as they are linearly stable in classical problem when 
< \l < fiRouth = 0.0385201. Further the n ormalizations of Hamiltonian and nonli near stability of £4(5) 
in the present of P-R drag has been studied by 



In this paper we generalized our previous paper 



Kushvah. Sharma. and Ishwari (l2007albllcl) 



Kushva h (2008a) in which linear stability has been 



examined in the generalized photogravitaional Chermnykh's problem, we have seen that the collinear 
points are linearly unstable and triangular points are stable in the sense of Lyapunov stability provided 
A* < URouth = 0.0385201. In present paper we have obtained the equations of motion and equilibrium 
points, zero velocity curves and the linear stability in the generalized photogravitaional Chermnykh's prob- 
lem. We have found that the collinear points deviate from the axis joining the two primaries, while the 
triangular points are not symmetrical due to Poynting-Robertson effect. We have examined the effect of 
gravitational potential from the belt, oblateness effect and radiation effect on the Lyapunov stability. 



2 EQUATIONS OF MOTION AND POSITION OF EQUILIBRIUM POINTS 



Let us consider the model proposed by lMiyamoto and Nagail (119751) . according to this model the potential 
of belt is given by: 

h 2 M b [ar 2 + (a + 37V)] (a + N) 2 



V(r,z) 



N 3 



r 2 + (a + N)' 



5/2 



(2) 



where M b is the total mass of the belt and r 2 = x 2 + y 2 , a, b are parameters which determine the density 
profile of the belt, if a = b = then the potential equals to the one by a point mass. The parameter a 
"flatness parameter"and b "core parameter", where N — \/ z 2 + b 2 , T = a + b, z 



V(r,0) 



M h 



Vr 2 + T 2 



0. Then we obtained 
(3) 



M bV Ac j n 

(r 2 +T 2)3/2 ■ AS in 



Kushva h (2008a bj), we consider the barycentric rotating 



and V x = (rB ffi )3/a , V y 

co-ordinate system Oxyz relative to inertial system with angular velocity u> and common z-axis. We have 
taken line joining the primaries as x-axis. Let m% , be the masses of bigger primary(Sun) and smaller 
primary(Earth) respectively. Let Ox, Oy in the equatorial plane of smaller primary and Oz coinciding with 
the polar axis of iri2- Let r e , r p be the equatorial and polar radii of mi respectively, r be the distance 
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of the masses and distance between primaries as unity, the unit of time i.e. time period of mi about m 2 
consists of 2ir units such that the Gaussian constant of gravitational k 2 = 1. Then perturbed mean motion 
n of the primaries is given by n 2 = 1 + + 2 2A ^ b 2 rc 3/2 , where r 2 = (1 — fJ,)q^ 3 + fx 2 ■ For simplicity, we 



set r 



2 1 ( r 2 +T 2)3/2, 

0.9999, T = 0.01 for further numerical results, where A 2 



5r 2 



is oblateness coefficient 



of m 2 . where u = — ^J 3 — is mass parameter, 1 — u = — ^ — with mi > mo. Then coordinates of mi 

^ ' mi+m2 * ' r° mi _)_ m2 ± z ± 

and m2 are (xi,0) = (—fi,0) and (£2,0) = (1 — /i, 0) respectively. Further, in our consideration, the 
velocity of light needs to be dimensionless, too, so consider the dimensionless velocity of light as Cd = c 
which depends on the physical masses of the two primaries and the distance between them. The mass of 
Sun mi w 1.989 x 10 30 kg w 332, 946m2(The mass of Earth), hence mass parameter for this system is 
H = 3.00348 x 10~ 6 . In the above mentioned reference s ystem the we determined the equations of motion 



of the in finitesimal mass particle in xy-plane. Now using 
( 2008afb ). then the equations of motion are given by: 



Miyamoto and Nagail d 19751) profile and 



Kushvah 



x — 2ny = U x — V x = Q x , 



y + 2nx = U y — V y 



(4) 
(5) 



where 



(1 - n)q\{x + fi) fJ,(x + fJ,-l) 



3 fiA 2 (x + fj, — 1) 
2 rf 



M b x 



r \ ( r 2 +T 2)3/2 
1 J { X + fi) ^ + ^ + yy]+x- 



Cl y = n y 



3 \iA 2 y 



(1 - fi)qiy fiy 



M b y 
(r 2 +T 2 ) 3/2 



ny 



^\^[(x + fi)x + yy}+y + n(x + fi) 



n 



n 2 {x 2 + y 2 ) (l-fx)qi fi 



n 



r-2 



M h 



^ + 

2r\ ( r 2 +T 2)V2 

(x + fi)x + yy 



2r{ 



n arctan 



V 



x + (1 



Wi 



(i-M)(i-gi) 



(6) 



qi = 1 — is a mass reduction factor expressed in terms of the particle radius 



a, density p radiation pressure efficiency factor \ ( m C.G.S. system):qi = 1 



5.6x10' 



-X- The energy 



integral of the problem is given by C = 2S1 — x 2 — y 2 , where the quantity C is the Jacobi's constant. 
The zero velocity curves C = 2£l(x, y) are presented in various frames of figure (Q~|i for the entire range 
of parameters 51,^2, M b details of the frames are given in the table ( [TJ. We have seen that in frame 
B(qi — .5) and C(qi = 1) there are closed curves around the £4(5) so they are stable but the stability 
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Table 1 Zero velocity curves when T = 0.01, fi = 0.025 



Frame 
when Mb = 0.2 


I(A 2 = 0.00 ) 
oval around L 4 fc\ 


H( A 2 = 0.02) 
oval around L\(k\ 


IU(A 2 = 0.04) 
oval around L\ik\ 


About the stability of 


A (qi = ) 
B ( qi = 0.5) 
C (gi = 1 ) 


No 

very small 
yes 


No 
very small 
yes 


No 
very small 
yes 


Unstable 
Stability reduced 
Stable 


Frame 
qi = 1A 2 = 0.02 


I(M 6 = 0.25) 


H(M b = 0.5) 


\M h = 0.75) 


Stability affected by 
mass of the belt 


D 


low effect of belt 


midium effect of belt 


very high effect of belt 


yes 



around Lm§\ disappeared so they are unstable due to P-R effect. When the P-R effect is absent as in frames 
D(q l = 1,A 2 = 0.02/ - M b = 0.25, II - M b = 0.5, III - M b = 0.75) the effect of oblateness and mass 
of the belt is presented, we have seen that the position of ovals around the £4(5) are different but still £4(5) 
are stable. 

The position equilibrium points are given by putting Q x — fl y = i.e., 

(1 - ju)<?i (x + n) fi(x + fi—1) 



n 2 x - 



3 [iA 2 (x + /j — 1) 

~2 4 

Wrny _ 
' o u, 



M b x 



(r 2 + T 2 f' 2 



2 (i-^)ais/ 
n y 3 — 



fiy 3 fiA 2 y 



2 r 2 5 



Win(x + jtt) 



(r 2 + T 2 ) 3/2 

when (Wi 7^ 0), from equations (UJ [8]l we obtained: 



1/3 

n = qi 



1 - 



nWi 



= 



2 



r2 



6(1 - n)y 

(1 - 2r c )M» (l-jjgfr 

3(r 2 +T 2 ) 3/2 

1 | M (l - 2r e )M & [ nW 1 
3 (r 2 + T 2 ) 3/2 



From above, we obtained: 



(1 - 2r c )Af b ( 



2/3 



1 



1 + 

2(1-/.) 



nWi 



2(1 

2/3 



+ 



3A 2 



± 



(r c 2 + T 2 ) 3/2 
nWi 



1/2 



. 5 - 2r )M b 

W 2 (r 2 + T 2 ) 3/2 



-2/3 



1/2 



(7) 



(8) 



(9) 
(10) 



(11) 



(12) 



From equations ( [7j [8]) the value of y is allwayes positive, hence the the equilibrium points are no-longer 
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then from equations (|4|i and (|5]l we obtained the triangular equilibrium points as: 



2/3 nWi 
2 6/*( 1 - n)y 



T ~ 2(1 - 14 



(l-2r e )M 6 [{l- ^}ql 
3 (rl + T 2 f /2 



2/3 



(13) 



y = ± 



2/3 
?1 



4 - g 2/3 + 2 (q 2/3 -2)4, 



2nl^i (q 2/3 - 2^ 
3yu(l - /i)y 



4(2r c - l)M b 



( 2/3 



-3 - 



3M 2 (^ /3 -3) 
2(1-M) 



1 1/2 



3 (r 2 + T 2 f /2 



(14) 



All these results are similar with 
others. 



Szebehelyld 19671) . 



Ragos and Zafiropoulosl(ll995l) . 



Kushvah(2008 



J3) 



and 



3 LYAPUNOV STABILITY 

In this section we will examine the P-R effect on the linear stability conditions. In mathematics, the no- 
tion of Lyapunov s stability occurs in the study of dynamical systems. In simple terms, if all solutions of 
the dynamical system that start out near an equ equilibrium m point L L stay near Li forever, then Li is 
Lyapunov's table. Let the position of any equilibrium point is (x*, y*) the taking x = x * +a, y = y* +(3, 
where a = £e At , (3 = rje xt are the small displacements £, 77, A are parameters, then the equations of 
perturbed motion corresponding to the system of equations (|U), ([5} are as follows: 

&-2n/3 = ani x +pn* xy + an* xi +l3n* xi (15) 
(3 + 2na = afi* + /?0* + afl* . + $(1*. (16) 



where superfix * is corresponding to the equilibrium points. 



(A 2 - xn* xx - n* xx )t + [-(2« + n* xy )x - n* xy ] v = o 
[(2n - n; ± )x - n; x ]t + (x 2 - xn; y - n;> = o 



(17) 
(18) 



this system has singular solution if, 



x 2 -xa xx -n* x -(2n + n* xi )x-n* 
(2n - n* yx )x - n* yx x 2 - xn; y - n; y 
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at the equilibrium points: 



a = 0^5—, 

r i* 

b = 2n 2 - /* - 
c = -a(l + e), 



3Af fc T 2 



-A, 



3M 6 



3/iA 2 

r i* T (r 2 +T 2 ) 5/2 

l* r 2* V 



5A 2 \ 2 
2r 2 „J y * 



(r 2 +T 2)5/2 



2Wf 



d = (n 2 - f,) 
+9/i(l - /i)y 2 
3M 6 



+ 2/* - -55 
9i 



3M 6 T 2 



1*' 2* 

,5 



r t (r 2 +T 2 ) 5/2 

(i-m)( 



( r 2 +T 2)5/2 I r ^ ^ 

6/inWiy* f (a;* + + n - 1) + y 2 



2M b (a;* + A*) + vl 
{rl+T 2 f 2 



where /* 



(l-M)gi 



3 A 2 
2 



( r 2_j i y f 2'-)5/2 ■ The points Li, L2, £3 no longer lie along the line 



joining the primaries, since the condition is not satisfied for them, so taking y — > 0, — — > because 

1/0 

y >> Wi,x >> W\, from ([9} we have ri « [■%] . Since f* > 1 and characteristic equation ( |T9b 
positive root for collinear points so they are unstable. Using Ferrari's theorem the roots of characteristic 
equation ([T5]) are given by: 



a + A 



- B 



(20) 



where A = ±V8Z - 46 + a 2 , and B = (§ + c*ia 2 ) (l ± yi + 8ai) T V1 1 + 8 e ai , c*i 
from this the characteristic roots are given by: 

a (l + VI + Sol) 
^1,2 — : 



(l+e)(l+e 2 -fc)+d 
2(fc 2 -4d) 



> 



^ w ^ 2 (i + VT+toT) 



^3,4 



16 

a (1 - VTTSoT) 



(21) 



±1 



(1- Vl + 8ai) 
16 



where 



Bo 



1 + e 



(22) 
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from equations( |2T1 122b . we found that at least one of the roots \(i — 1, 2, 3, 4) have a positive real part 
due to P-R ef fect thus the triang ular equilibrium point s are unstable in the sense of Lyapunov stability result 



are similar to 



Chernikov(1970) and 



Kushvah (2008b). 



4 CONCLUSION 

We have seen the collinear points L%, L 2 , £3 no longer lie along the line joining the primaries and they are 
unstable in classical case. If q± = .5, q% = 1 there are closed curves around the £4(5) so they are stable but 
the stability range reduced[see frame B(qi = 0.5) in figure ([TJ] due to P-R effect. When the P-R effect is 
absent as in frames D(q 1 = 1,A 2 = 0.02/ - M b = 0.25, II - M b = 0.5, III - M b = 0.75) the effect of 
oblateness and mass of the belt is presented, we have seen that the positions of ovals around the £4(5) are 
different but still £4(5) are stable. In frame A(qi = 0) the closed curves around £4(5) disappeared so they 
are unstable due to P-R effect. 
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Fig. 1 This figure show the zero velocity curves for /i = 0.025, T = 0.01, frames (A-I to 
III)qi = 0.00, M b = 0.02, A 2 = 0.00 - -0.04, (B-I to III) 9 i = 0.50, M b = 0.02, A 2 = 
0.00 --0.04,(C-I to III) g i = 1.00, M b = 0.02, A 2 = 0.00- -0.04,(D-Ito lll)qi = 1.00, M b = 
0.25,0.5,0.75,^2 = 0.02 



